Compact Boson Stars 



Betti Hartmann", Burkhard Kleihaus 6 , Jutta Kunz 6 and Isabell Schaffer" 

"School of Engineering and Science, Jacobs University, Postfach 750 561 

D-28725 Bremen 
b Institut fur Physik, Universitat Oldenburg, Postfach 2503 
D-26111 Oldenburg, Germany 

March 3, 2013 

cn" 

i— H 

o 

^~». Abstract 

We consider compact boson stars that arise for a V-shaped scalar field potential. They 
represent a one parameter family of solutions of the scaled Einstein-signum-Gordon equations. 
We analyze the physical properties of these solutions and determine their domain of existence. 
Along their physically relevant branch emerging from the compact Q-ball solution, their mass 
increases with increasing radius. Empoying arguments from catastrophe theory we argue that 
this branch is stable, until the maximal value of the mass is reached. There the mass and size 
are on the order of magnitude of the Schwarzschild limit, and thus the spiralling respectively 
»' j . oscillating behaviour, well-known for compact stars, sets in. 

1 Introduction 

> . 

, When a complex scalar field is coupled to gravity, boson stars arise as stationary localized solutions 

Q\ • of the coupled Einstein-Klein-Gordon equations [TJ[21[3]. The physical properties of the boson stars 

depend strongly on the scalar field potential (see e.g. the review articles [H El El [5] ) . 

If only a mass term is present, for instance, but no self-interaction, so-called mini boson stars 
, arise. Their masses are relatively small, being bounded by a maximal mass M max on the order 

of the Planck mass squared, divided by the boson mass, Mp^/ms When a repulsive quartic self- 
interaction is included larger boson stars are obtained, having a maximal mass M max on the order 
of V\M^/m%, where A is the self-coupling constant [9^ . The presence of a sextic potential, on the 
other hand, leads to boson stars with a maximal mass M max on the order of Mpjrrv^ [4]. These 
solitonic boson stars even possess a flat space limit, where they correspond to non-topological 
solitons [TO: (or Q-balls [TT]). 

All these types of boson stars, whether small or large, do not possess a sharp radius. The 
scalar field exhibits an exponential fall-off as dictated by the particle mass, but this does not yield 
a unique boundary for the star. Thus a large variety of radius definitions are employed, when one 
is interested in extracting the relation between the mass of the boson stars and their size. 

Here we consider boson stars, whose radius is uniquely defined, because these solutions are 
compact. Their scalar field vanishes identically outside the radius of the star. These compact 
boson stars arise when a V-shaped self-interaction potential is employed. In the flat space limit 
this potential allows for compact Q-balls, which represent solutions of the signum-Gordan equation 
for the scalar field [T2j [13] . Interestingly, when the electromagnetic field is coupled, even shell- like 
solutions appear, where the scalar field vanishes identically outside a finite shell [TBI UM [T5] , 

By solving the coupled Einstein-signum-Gordon equations, we determine the domain of exis- 
tence of these compact boson stars. In particular, we show that by scaling the equations appropri- 
ately, the solutions depend only on a single parameter, varying between zero and a finite maximal 
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value. We then analyze the physical properties of these compact boson stars and address their 
stability from a catastrophe theory point of view [16j [TT1 [TBI [19l [20] . 

The paper is organized as follows. In section 2 we present the action, the Ansatz, the equations 
of motion together with the scaling property, and the global charges. We present our solutions 
in section 3, and discuss their physical properties. We end with a conclusion and an outlook in 
section 4. 



2 Action 



We consider the action of a self-interacting complex scalar field $ coupled to Einstein gravity 

- 4 h 2 



s 



r — (d^y (a***) - t/(|$o 



16nG mo 



-gd 4 x, 



(1) 



with curvature scalar R, Newton's constant G, a mass scale mo, and the asterisk denotes complex 
conjugation. The scalar potential U is chosen as 



U{\®\) = 2A|$| 
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The dimensions of $ and A are 1/length 3 / 2 , respectively [energy] /length 3 / 2 . Next we introduce 
the scaled scalar field \ff = fij Jm§Q. The resulting action reads 
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where we defined A = ^/m Q /h\. 

Variation of the action with respect to the metric and the matter fields leads, respectively, to 
the Einstein equations 
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with stress-energy tensor 
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and the matter field equation, 
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where V M denotes the covariant derivative. 

To construct static spherically symmetric solutions we employ Schwarzschild-like coordinates 
and adopt the spherically symmetric metric 



with 



ds 2 = g^dx^dx" = -A 2 Nc 2 dt 2 + N^dr 2 + r 2 (d6 2 + sin 2 6d(j) 2 ), 
N = 1 2Gm(r)^ 



The Ansatz for the matter fields has the form 
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The resulting Einstein and field equations are 
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In the next step we introduce the dimensionless coordinate x = r/oj and the dimensionless scalar 
field h = u?IM>. This yields 



V! = x*\Nh" + -^h' + 2h) , (13) 
A' = 2ax(Ah' 2 + ^h^ , (14) 

h " = -^ h+ ^W-^ff- h '> (15) 



where 



a = — ^, (16) 

cj 5 c 2 ujG , , 

are the dimensionless coupling parameter, respectively mass function and N = 1 — 2oyx/x. 

Let us now specify the boundary conditions for the metric and matter functions. For the metric 
function A we adopt 

A(x ) = 1 , (18) 

where xq is the outer radius, thus fixing the time coordinate. For the metric function N(x) we 
require 

N(0) = 1 , (19) 
and for the scalar function we require at the origin and at the outer radius xq, respectively, 

h'(0) = , h(xo) = , h'(x ) = . (20) 

Note that h(x) = 0, A(x) = l,[i(x) = const is an exact solution of the equations (|13 p -(|15 p in 
the exterior region x > xq. 

Hence the mass of the boson star is given by 

2 2 

M = ™M(oo) = -%l<*o) . (21) 
uj G U) G 

and the areal radius of the boson star is 

R = x q /uj . (22) 

As a consequence, ^§|jr = 2a . 

The particle number Q is computed from the density 

lh fl " . (23) 
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Using the ansatz and dimenensionless quantities yields 



P = 



A 2 1 h 2 
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Integration then gives the particle number 




py/^gdrddd(j> 



A^4tt 

w 6 he 




(25) 



where 
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For the central energy density we find 



^ 2 (hi „, 



) 




(27) 



where ho = h(0), Aq = A(0) and the boundary condition h'(0) — were used. The central mass 
density is then defined as no = £o/c 2 - 

3 Results 

We solve the coupled system of equations numerically employing a Newton-Raphson scheme. 

Let us begin our discussion of the results with the flat space-time limit, the compact Q-ball. 
Here the scalar field function can be obtained analytically [TH [T3] 



where xq satisfies xq — tanxo, i.e., xq ~ 4.4934, while the metric functions are trivial, N(x) = 
A(x) — 1. This analytical solution is our starting point, obtained in the limit a — > 0. It is seen in 
Figs, la-lc. 

As we increase a from zero we couple gravity, since a ~ G. With increasing a a first branch 
of compact boson star solutions evolves, where the scalar field and the metric functions change 
continuously. This first branch ends at a finite maximal value of a, a max - This is seen in Fig. Id, 
where we exhibit the value of the scalar function at the origin, /i(0) versus a. At a max a second 
branch bends backwards towards smaller values of a. While h(0) now increases monotonically, a 
starts to exhibit damped oscillations, tending towards a limiting value, au m . 

The scalar field and metric functions of this family of solutions are exhibited in Figs, la-lc for 
a representative set of values of a, as marked by the dots in Fig. Id. In the oscillating regime of a 
the functions start to show a distinct behaviour, deviating from the smooth first branch pattern. 
In particular, the function h(x) increases steeply at the center of the boson star, while the function 
A(x) at the same time decreases strongly at the center, tending towards zero in the limit a — > an m . 

Let us next consider the relation between the parameter a and the boson star radius parameter 
xq. We exhibit a versus xq in Fig. 2a. The first branch starts from the flat space-time limit, a = 0, 
where a increases with decreasing xq . The previously seen oscillating behaviour of a for the higher 
branches now translates into a spiralling behaviour, when considered versus the radius parameter 

X . 

Likewise, the mass parameter /x(oo), presented in Fig. 2b, exhibits an analogous spiralling 
behaviour, except that along the first branch the mass parameter is decreasing with decreasing xq, 
and the spiral is inverted. Clearly, a limiting value /in m of the mass parameter is approached at the 
center of the spiral, where a — > an m . The maximal value of /i corresponds to the flat space-time 
limit, marked by a dot in the figure. 
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Figure 1: The scalar field function h(x) (a), the metric function N(x) (b), and the metric function 
A(x) (c) versus x for several values of a. Also shown are the values of the scalar field function at 
the origin, h(Q), versus a (d). Here the dots mark the set of values of a selected in (a)-(c). 



In Fig. 2c we exhibit the mass parameter /i(oo) versus the dimensionless central density e c . 
Here the mass parameter exhibits damped oscillations, as it approaches /in m . Finally, in Fig. 2d, 
we consider the dependence of the mass parameter on the scaled particle number. Along the first 
branch, starting from the flat space-time limit, the relation between the two quantities is almost 
linear. Then again a spiralling behaviour sets in. 

These figures with the dimensionless quantities represent the full domain of existence of the 
compact boson star solutions. All physical solutions with dimensionful quantities can be obtained 
from these dimensionless solutions by appropriate scaling. 

To make contact with real stars and astrophysics, let us now discuss such sets of physical 
solutions. To this end, we now consider the solutions for fixed models, i.e., for fixed values of the 
coupling constant A. Moreover, we take the physical value for Newton's constant. Then the value 
of a yields the value of the frequency u). Together with /^(oo) we then obtain the corresponding 
value of the mass M of the respective compact boson stars. Likewise, by employing the radius 
parameter xq together with the frequency u> we obtain the physical value for the areal radius R 
of the respective compact boson stars. 

We exhibit the mass M in units of the solar mass Mq versus the radius R in kilometers in 
Fig. 3a for several values of the coupling constant A, determining the strength of the potential, 
i.e., for several fixed models. We observe, that the maximal value of the mass and the maximal 
value of the radius increase with decreasing coupling strength A of the potential. 
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Figure 2: Coupling a versus radius parameter xq (a); mass parameter [i versus radius parameter 
xq (b), central density e c (c) and scaled particle number Q/Qo (d). The dots mark the flat limit. 



Interestingly, the compact boson stars also correspond to physically compact stars, since for 
a mass on the order of the solar mass their radius is on the order of ten(s) of kilometers, thus 
corresponding in mass and size to neutron stars. However, their mass increases with increasing 
size, except in a small region close to the maximal value of the radius. This behaviour is in contrast 
to neutron stars, where we observe an increase of the mass with decreasing radius. Instead, the 
observed dependence of the mass on the radius is more in line with quark stars [3TJ [231 HS1 HI] • 
But such a behaviour has also been seen in boson stars of the soliton type [3D] . 

As seen in Fig. 3a, after reaching the maximal value of the mass the solutions show the 
well-known spiralling behaviour of highly compact stars. For neutron stars this behaviour is 
encountered, when the mass and size of the stars approach the values of the Schwarzschild black 
hole solution. Therefore, for comparison, we have also included the Schwarzschild limit in the 
figure. While we are in the vicinity of this limit, when the spirals are encountered for the compact 
boson stars, we do not get very close. In contrast, the soliton type boson stars practically reach 
the black hole limit [20] . 

It is also interesting to consider the dependence of the mass M on the central mass density no 
of the boson stars. Therefore we exhibit this dependence for the same set of solutions in Fig. 3b. 
First, the mass reaches monotonically its absolute maximum, and then exhibits the well-known 
oscillating behaviour, also observed in neutron stars. However, we note that the higher the mass 
at the maximum, the lower the central mass density of the compact boson stars. We attribute 
this behaviour to the necessity of lowering the potential strength in order to reach higher masses. 
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Figure 3: The mass in units of the solar mass versus the areal radius (a) and the central mass 
density (b) for several values of the potential strength A (Aq = y/2 x 10 30 kg m~ 3 s~ 2 ) . 



Let us finally address the stability of the configurations. A mode analysis, as performed in the 
flat space limit [55] is very involved. Therefore we will use arguments from catastrophe theory to 
discuss the stability [TBI HZ1 HH HH1 12D] ■ We will base our discussion on Fig. 3b, where the mass is 
shown versus the central mass density for several values of the potential strength. In particular, 
we choose the mass and the potential strength as the two control parameters, and we choose the 
central mass density as the single behavior variable [HI \19\ [20] . 

To analyze the stability of the compact boson stars, we start from the equilibrium space M. = 
{no, M, A}, as exhibited in Fig. 3b. According to catastrophe theory, the stability changes only at 
the turning points, where dM /dug = (for fixed values of A). Thus passing a turning point means 
changing the stability of the boson star configurations. The first branch starts from a stable [21)] 
solution in flat space-time. Thus this branch should be stable, until the maximal value of the mass 
is encountered. At the maximum of the mass, however, stability should change, and therefore the 
solutions in the spiral should be unstable. In fact, we expect that a mode analysis will reveal, 
that the solutions will become increasingly unstable inside the spiral, collecting collecting more 
and more unstable modes. 



4 Conclusions and Outlook 

We have studied compact boson stars obtained from a V-shaped interaction potential. Thus the 
scalar field is confined to a finite region. By analyzing the scaled set of field equations, we have 
obtained the full set of solutions, consisting of two parts. The first part corresponds to the branch 
of solutions that emerges from the flat space-time solution, and the second part, depending on the 
physical quantities studied, consists of a spiralling or oscillating set of solutions. 

By specifying the constants of the model, we have translated these dimensionless results into 
families of compact boson star configurations whose masses and sizes can be compared to real 
astrophysical objects. Along the physical branch, emerging from the flat space-time solution, the 
mass increases with the radius, until a maximal value of the radius is reached. This behaviour is 
unlike neutron stars and more akin to quark stars. But as for these fermionic stars, the size of 
solar mass compact boson stars in on the order of ten(s) of kilometers. 

However, the masses of these compact boson stars can be made arbitrarily high, and their sizes 
increase accordingly, when the potential strength is made sufficiently small. These compact boson 
stars thus could represent huge compact astrophysical objects, not too far from their corresponding 
Schwarzschild black hole limit. 

Employing arguments from catastrophe theory, we have argued that the physical branch of 
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these solutions, that emerges from the stable flat-space solution, is stable as well, while the re- 
maining configurations are unstable. 

As a next step we will include rotation [28 . Rotating boson stars have been obtained before 
for non-compact stars [371 IHJ HH 1301 I3H 132] ■ m particular, they exhibit a quantization relation 
between the angular momentum and the particle number. We expect, that the inclusion of rotation 
in these compact boson stars might lead to interesting new observations. 

It should also be interesting to construct interacting compact Q-balls and boson stars, which 
should arise in the presence of several complex scalar fields [321 1331 131] • 
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